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Abstract 

We study the roughening transition of the dual of the 2D XY model, of the 
Discrete Gaussian model, of the Absolute Value Solid-On-Solid model and of the 
interface in an Ising model on a 3D simple cubic lattice. The investigation relies on a 
renormalization group finite size scaling method that was proposed and successfully 
tested a few years ago. The basic idea is to match the renormalization group 
flow of the interface observables with that of the exactly solvable BCSOS model. 
Our estimates for the critical couplings are /3jf y = 1.1199(1), K^ G = 0.6653(2) 
and Kft SOS = 0.80608(2) for the Ay-model, the Discrete Gaussian model and 
the Absolute Value Solid-On-Solid model, respectively. For the inverse roughening 
temperature of the Ising interface we find K^ mg = 0.40758(1). To the best of our 
knowledge, these are the most precise estimates for these parameters published so 
far. 



1 Introduction 



Among the phase transitions that occur in 2D or effectively 2D statistical systems, those of 
the so called Kosterlitz-Thouless (KT) type belong to the most challenging. The KT phase 
transition is of infinite order: The free energy and all its derivatives stay finite at the transition 
point. Despite the relatively simple arguments that suggest the existence of such a transition in 
a variety of systems, a rigorous proof of the KT nature of the phase transition in many physically 
interesting systems is still lacking. Also most of the numerical studies (many of them Monte 
Carlo studies of the 2D XY model), could not provide an unambigious confirmation of the KT 
scenario. A number of references will be given in section || 

The reason for the problem is the appearance of corrections to scaling that vanish only 
logarithmically with the system size. Most of the investigations based on simulations of KT 
models on finite lattices suffer from these corrections. 

Note, however, that there exists at least one 2D lattice model, which has been proven 
to undergo a KT transition by exact solution. This is the Body Centered Solid-On-Solid 
(BCSOS) model [|J, the configurations of which are, up to boundary conditions, in one-to-one 
correspondence to those of a special 6- vertex model ||, [|, ||], the F-model. 

A few years ago we proposed a method that allows to investigate the KT transition of a 
given model by comparing its block spin renormalization group (RG) flow (on finite lattices) 
with that of the BCSOS model |^, [?]]. A matching of the two RG flows at long distance (large 
blocks) demonstrates that the models belong to the same universality class, which is the KT 
class here. 

In contrast to the usual approaches, the method introduces systematic errors that decay 
like L~ 2 , where L is the size of the lattices involved in the computations. 

Our approach has been successfully applied to the Absolute Value Solid-On-Solid (ASOS) 
model, the Discrete Gaussian (DG) model and the dual of the standard XY model in two 
dimensions §. Successful applications to the interface of the Ising model were performed 
in M, and, recently, in H. 

In the present paper we improve on the results of ref. || by using larger lattice sizes 
and increasing the statistics by a factor of about 100. This became affordable both by the 
availability of faster computers and the use of more efficient program code. 

This paper is organized as follows: In section 2 we define the models and state the exact 
results for the BCSOS model relevant to our study. We briefly discuss the KT flow equations. 
Section 3 is devoted to a description of the matching method. In section 4 we present and 
discuss our numerical results. A comparison with previous estimates of the critical couplings 
and non-universal parameters is presented in section 5. Conclusions and an outlook follow. 

In the preprint version of this paper, part of the tables is presented in a section after the 
Bibliography. 

2 Ising Model Interfaces and Solid-On-Solid Models 
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2.1 Ising Model Interfaces 

We consider the 3D Ising model on the simple cubic lattice, with Hamiltonian 

H = - S x S V , S x = ±1 . (1) 

<x,y> 

The sites of the lattice are labelled by integer coordinates x = (xi,X2,x^). The sum in eq. (|l]) 
is over all (unordered) nearest neighbour pairs of sites in the lattice. The partition function is 

Z = £ exp (-K'H) . (2) 

M 

Here, the summation is over all possible configurations of the Ising spins. The pair interaction 
is normalized such that K 1 = l/(fc#T), where ks denotes Boltzmann's constant, and T is the 
temperature. 

At a critical coupling K\ (the estimate of a recent study ||] is = 0.2216546(10)) the 
infinite volume limit of the model undergoes a second order phase transition. For K > K^, 
the system shows spontaneous breaking of the reflection symmetry. 

In order to study interfaces separating extended domains of different magnetization, we 
consider lattices with extension L in the x\- and ^-directions and with extension D in the 
^-direction. We generalize eq. ([[]) to 

H = S ' k X y S X Sy . (3) 

<x,y> 

The lattice becomes a torus by regarding the opposite boundary planes as neighbour planes. 
For the Ising spins s we will apply antiperiodic boundary conditions in X3-direction, by letting 
k X y = — 1 for the links that connect the uppermost with the lowermost plane. For the other 
links we set k xy = 1. 

For sufficiently large K 1 and large enough L, the imposure of antiperiodic boundary con- 
ditions forces the system to develop exactly one interface, a region where the magnetization 
rapidly changes sign. This interface is parallel to a (001) lattice plane. 

The Ising (001) interface undergoes a roughening transition at an inverse temperature = 
l/(&B^i?) that is nearly twice as large as the bulk transition coupling K l c given above.[] In this 
work, we shall determine a new estimate for K^, and also for other parameters of the roughening 
transition. For a collection of previous estimates, see section ||. 

At the roughening transition, the large scale interface behavior changes from being rigid or 
smooth at low temperature to being rough at high temperature. The transition shows up in a 
characteristic behaviour of various quantities. For example, in the smooth phase, the interfacial 
width stays finite when L tends to infinity, while it diverges logarithmically with the system 
size in the rough phase [11, [l2j ]. For general introductions to roughening, see [jl^, 14, 15|. For 



comparisons of real life experiments with theory see, e.g., refs. [16 



1 A pioneering work on this issue is ref. uXM 
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2.2 Solid-On-Solid Models 



A fairly good approximation of the Ising interface is given by Solid-On-Solid (SOS) models to 
be introduced in this section. The SOS approximation amounts to ignoring overhangs of the 
Ising interface and bubbles in the two phases separated by the interface. For a review of exact 
results on SOS type of models, see, e.g. By duality [17] and other exact transformations 



(see, e.g., jig] ), SOS models have been shown to be equivalent to a variety of other statistical 
models. 

All SOS models that we shall consider have in common that they are 2D lattice spin models. 

Our first example of an SOS model is the Absolute- Value-Solid-On-Solid (ASOS) model. 
It can be considered as the SOS approximation of an (001) lattice plane interface of an Ising 
model on a simple cubic lattice. The model is defined by the Hamiltonian 

Hasos = K ASOS \h x - h y \ . (4) 

The spin variables h x take integer values. Here and in the following, the Boltzmannian will 
always be exp(— H). A factor 1/(&bT), where &b denotes Boltzmann's constant and T the 
temperature, is absorbed in the definition of the Hamiltonian. 

We interprete the h x as heights with respect to a certain base. For finite positive K ASOS 
the Hamiltonian will favour that neighbouring spins take similar values. When K ASOS is large 
enough, the surface will not fluctuate too wildly (smooth phase). On the other hand, if K ASOS 
is below a certain critical value, the surface becomes "rough", and, e.g., the surface thickness 
diverges when the system size goes to infinity. 

Let us now turn to the Discrete Gaussian (DG) model. The Hamiltonian is 

H DG = K DG (h x - h y f . (5) 

The spin variables h x take integer values. Note that the Hamiltonian looks exactly like that of 
a continuous Gaussian model. However, the restriction of the h x to integer values introduces 
a nontrivial interaction. The Discrete Gaussian model is dual to the XY model with Villain 
action u7l.n This model is defined by the partition function 



z v = n de - n B(e x -o y ), (6) 

"'~ 7r x <x,y> 

with 

oo 

B(0) = Yl exp(-I/3 y (0-2^) 2 ) (7) 

p=— oo 

and 

' K DG . (8) 



Wv 



What is called Hamiltonian in the language of Statistical Mechanics is called action in the framework of 
Euclidean Quantum Field Theory. 



3 



The index "V" here refers to "Villain". 

The XY model with "standard (cosine) action" has the partition function 



z X Y = r n d0 - ex p [p xY e cos ( q - - °y) ) 

J -' K x \ <x,y> J 



0) 



The standard action is the mostly discussed action for an XY model. The dual of this model 
is given by the partition function 

z s x°y s =e n 'iw^ y )> (io) 

{h} <x,y> 

where the I n are modified Bessel functions. Again h x is integer. 

We finally introduce the Body Centered Solid-On-Solid (BCSOS) or F-model. The BCSOS 
model was introduced by van Beijeren [|j as an SOS approximation of an interface in an Ising 
model on a body centered cubic lattice on a (001) lattice plane. For a detailed analysis of 



this model with respect to roughening and surface structure, see [14, 19, 15]. The effective 
2D lattice splits in two sublattices like a checker board. In the original formulation, on one 
of the sublattices the spins take integer values, whereas the spins on the other sublattice take 
half-integer values. We adopt a different convention: spins on "odd" lattice sites take values of 
the form 2n + i, and spins on "even" sites are of the form In — |, n integer. The Hamiltonian 
of the BCSOS model can be expressed as 

H B csos = K BCSOS Y / \hx-h y \. (11) 

[x,y] 

The sum is over next-to-nearest neighbour pairs [x,y], and nearest neighbour spins h x and h y 
obey the constraint \h x — h y \ = 1. Van Beijeren [Q] showed that the BCSOS model can be 
transformed into the F-model, which is a special six vertex model. The F-model can be solved 
exactly with transfer matrix methods H, |], 0]. The roughening transition occurs at 

K BCSOS = l ln2 (12) 



For K \ .fTft, the correlation length behaves like 




: bcsos „ i _ I _^ 1 K= (13) 



2.3 Renormalization Group Flow of Interface Models 

It is believed (though not proven rigorously) that, in the vicinity of the fixed point relevant for 
the KT transition, the RG flow of SOS models and also of the 3D Ising model interface is well 
described by two parameters f3 and z The two parameters are the inverse temperature and 
a fugacity z. 
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The 2D Sine Gordon model is especially suited to discuss the flow of these parameters with 
the length scale, since this model contains (3 and z as bare parameters in its Hamiltonian: 



rSG = Y& 2 - - ^E cos ( 2 ^) . ( 14 ) 

* <x,y> x 

where the <f> x are real numbers. For the continuum version of the model, with a momentum 
cutoff, one can derive the parameter flow under infinitesimal RG transformations jjj . It is given 
by 

x = — z 2 , z = — xz , (15) 

where z = const ■ z and x = 7r/3 — 2. const depends on the particular cutoff scheme used. The 
derivative is taken with respect to the logarithm of the cutoff scale. 

For large x the fugacity z flows towards z = 0. The large distance behaviour of the model is 
therefore that of a Gaussian model (without a mass term). For small x, z grows with increasing 
length scale. The theory is therefore massive, i.e., has finite correlation length. The critical 
trajectory separates these two regions in the coupling constant space. It ends at a Gaussian 
fixed point characterized by x = or (3 = ^. On the critical trajectory the fugacity vanishes as 

<t) = , (16) 

z +t 



where t is the logarithm of the cutoff scale. Eqs. fllq ) are the basis for KT theory. Its immediate 
consequences are derived in statistical mechanics text books, see e.g. |2(J. E.g., the correlation 
length in the smooth phase of an SOS model should diverge like 

£~ A exp (Ck- 1 ' 2 ) , «=*sfa, (17) 

when K — » Kr. We would like to emphasize another important consequence of the KT equa- 
tions that becomes apparent from the solution eq. (|l6|): At criticality, the fugacity, which 
parametrizes the deviation of the theory from a massless Gaussian model, decays with increas- 
ing scale t = InL, only like (InL) -1 . In lattice studies, L is more or less the lattice extension. 
Therefore, any method that is based on an observation of the Gaussian behaviour at long 
distance, suffers strongly from finite fugacity corrections even on very large lattices. 



3 The Matching Method 

The method of ref. || is closely related to the finite size scaling methods proposed by Nightin- 
gale |2l| and Binder |22| . No attempt is made to compute the RG flow of the couplings 
explicitly, but rather the RG flow is monitored by evaluating quantities that are primarily 
sensitive to the lowest frequency fluctuations on a finite lattice. One should stress that the 
method does not use any of the quantitative results of KT theory. Merely the qualitative result 
that there are two important coupling parameters in the flow is used. 
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In order to separate the low frequency modes of the field a block spin transformation [23, 24] 
is used. Blocked systems of size I x I are considered. The size B of a block (measured in units of 
the original lattice spacing) is then given by B = L/l, where L is the linear size of the original 
lattice. The linear blocking procedure defined by 

(t>x = B- 2 J2 h x , (18) 

where X labels square blocks of a linear extension B, is used. This linear blocking rule has the 
half-group property that the successive application of two transformations with a scale-factor 
of B have exactly the same effect as a single transformation with a scale- factor of B 2 . 

Motivated by the perturbation theory of the Sine Gordon model two types of observables 
are chosen: such that are "sensitive" to the flow of the kinetic term (flow of K), and such that 
are sensitive to the fugacity (periodic perturbation of a massless Gaussian model). For the first 
type of observables 

Ai = {{fix - 0y) 2 ) , (19) 
where X and Y are nearest neighbours on the block lattice, and 

A 2 = {(4>x - 4>zf) , (20) 

where X and Z are next to nearest neighbours, are chosen. Note that these quantities are only 
defined for I > 1. As a monitor for the fugacity the following set of quantities (defined for 
I = 1,2,4 and 8) is taken: 

A 3 = (cos(l • 2n4> x )) , 

A 4 = (cos(2 • 2vr^x)} • (21) 



3.1 Determination of the Roughening Coupling 

There are two parameters which have to be adjusted in order to match the RG flow of an 
SOS model or of the Ising interface with that of the critical BCSOS model: The coupling K 
of the Solid-On-Solid or Ising model and in addition the ratio b = B s /B B = L s /L B of the 
lattice sizes (and hence the block sizes) of the SOS or Ising model and the BCSOS model. In 
general a b ^ 1 is necessary to compensate for the different starting points of the two models 
on the critical RG-trajectory [|| . For the proper values of the roughening coupling and the 
matching constant b observables of the SOS and the BCSOS model match like 

Af tl (bB,K%) = A B (B,K%) + 0(B-n, (22) 

where i labels the observable and I the size of the blocked lattice. The 0{B~ UJ ) corrections are 
due to irrelevant operators, oj is the correction to scaling exponent. The perturbation theory 
of the Sine Gordon model suggests u> = 2. 

In order to obtain numerical estimates for the roughening coupling and the matching 
factor b for a given lattice size L B of the BCSOS model we require that eq. ( |2^ ) is exactly 
fulfilled for two block observables. 
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We solve the system of two equations for the two observables A^i and Ajj numerically by 
first computing the Kf t (b) and Kj^b) that solve the single equations for a given value of b. 
The intersection of the two curves Kf^b) and K^^b) gives us then the solution of the system 
of two equations. For an illustration of this method see figures 5 and 6 of ref. ||. 

In H we demonstrated, that the corrections to scaling for the observables A\ and A2 for 
SOS models are similar to those in the massless continuous Gaussian model. Therefore we 
considered the "improved" observable D\ which is defined as follows: 

' l(L) = ^ il(i) - (23) 
A^ is computed for the massless Gaussian model defined by 

Ho = \ E (^-VV) 2 - (24) 

<x,y> 

An improved quantity D2 is defined analogously. Explicit results for and are given 
in table ||. 

Obviously this modification does not affect the large L behaviour since 
0(L~ 2 ). It turns out that the results for our largest lattice sizes are virtually unaffected by 
this kind of improvement. 

3.2 Determination of Non-Universal Constants 

The matching programme also allows to determine the non-universal constants appearing in 
formulae describing the divergence of observables near the roughening transition. In Q we 
showed that the two non-universal parameters A and C determining the critical behaviour 
of the correlation length, cf. eq. (|l7|), can be determined from information of the matching 
procedure. For one of the models matched with the BCSOS model, one finds 

A sos = b m A BCSOS (25) 
c sos = q -i/2 C Bcsos (26) 

where the parameters A BCSOS and C BC can easily be extracted from eq. (|l3|). If 

BAfcsos /QAii 

dK BCSOS I Q K ' 

where quantities have to be taken at the roughening couplings, is the same for all observables, 
which is the case for our data, then q is given by 

K sos 

q = K^ sos R ' ^ 



For a more detailed discussion see ref. ||. 
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4 Numerical Results 



We simulated the BCSOS model at its critical coupling Kj^ csos = ^ln2 using the loop- 
algorithm of Evertz, Marcu and Lana |^| . One has to note that periodic boundary conditions of 
the F-model do not correspond to periodic boundary conditions of the BCSOS model. Therefore 
updates of loops that wind around the lattice are forbidden. 

We performed 10 7 measurements for all lattices sizes considered. We have chosen the 
number of loop-updates between two successive measurements such that the autocorrelation 
times were about 1. 

In addition to the observables An we measured the interface thickness, the total energy E 
and An x E, which is needed to compute derivatives of the observables with respect to the 
coupling. In order to save disc space we accumulated 1000 measurements before writing to the 
file. The statistical errors were computed by jacknifing the 10 prebinned data. As random 
number generator we used a combination of three shift register generators. 

We checked the reliability of the updating program by comparing the estimates from 10 8 
measurements for L = 4 with the exact results obtained by explicitly averaging over all BCSOS 
configurations. Our data are also consistent with those of ref. |J. In || lattices of size up to 
L = 96 were used, 4 x 10 6 measurements were performed, and the G05CAF random number 
generator of the NAGLIB was used. Note also that the computer programs of H were written 
independently of the programs used in the present study. 

The results for the BCSOS observables are summarized in tables |9] and 10. Our estimates 
for the slopes of the observables are given in tables 11 and ITTJ. With slope we here mean the 
derivative of the observables with respect to the coupling K BCSOS , taken at the critical value. 

We then performed the simulations for the ASOS, the DG and the dual of the XY model. 
The simulations of the ASOS and DG model were done using a demon version pq| of the 



Valleys-to-Mountains reflection (VMR) algorithm [27]. The simulation of the dual XY model 
was done using the standard version of the VMR algorithm. In both cases we used the G05CAF 
routine of the NAGLIB as random number generator. 

Again we performed 10 7 measurements and accumulated 1000 measurements before writing 
to the file. In order to get estimates for the observables in a neighbourhood of the simulation 
point, we employed a second order Taylor expansion (note that the prebinning forbids the 
use of a reweighting technique). We thus computed the first and second derivatives of the 
observables, 

A A 

— = (H)(A)-(HA), (29) 
(H 2 A) - 2(H) (HA) + 2(H) 2 (A) - (H 2 ) (A) , (30) 



and 

d 2 A 



dK 2 

for the ASOS and the DG model. In the XY case analoguos formulae were derived. We 
carefully checked by comparing with results obtained from simulations at shifted couplings 
that the Taylor expansion of the An to second order was sufficiently precise. 
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We performed the simulations at the previously best known estimates for the roughening 
couplings |], namely K^ G = 0.6645, K^ sos = 0.8061 and (3$ Y = 1.1197. 

In order to keep the length of this paper within reasonable bounds, we present the numerical 



results only for the XY model, see tables |13| and |14j. The tables for the other models are 
available from the authors upon request. 

Finally we performed the simulations for the Ising model. The simulations were done using 
the (VMR) algorithm [27[ adapted to the Ising interface as discussed in refs. [28, 12]. Following 



the proposal of ref. |29j we encoded all D spins with the same x± and %2 in a single word. We 
found no way to execute operations in parallel on the spins in one word. The advantage of the 
this particular coding merly lies in considerably reduced memory access time. 

The local update of the system was speeded up by replacing the standard Metropolis update 
by a demon update implemented in multispin coding technique. Also the elimination of the 
bubbles in the bulk phases was done using a multispin coding technique. 

In total these improvements led to a reduction of the CPU time required for a given statistic 
by a factor of about 4 compared with the code used in refs. |p8j , [|]. 

We performed the simulations at the previously best known estimate for the roughening 
couplings |g, Kl = 0.40754. 

We performed 3 x 10 6 to 8 x 10 6 measurements for lattice sizes ranging from 32 x 32 x 31 to 
192 x 192 x 31. Again we pre-binned the results of 1000 measurements before writing to disc. 

For the Ising interface we also computed the Taylor expansion of the observables to second 
order. 

In the case of the Ising interface we performed, in addition of the matching with the BCSOS 
model, the matching with the ASOS model at K^ sos = 0.80608, which is our present estimate 
of the roughening coupling for the ASOS model. For this purpose we performed additional 
simulations of the ASOS model at K^ sos = 0.80608 for the lattice sizes 24,48, 56, 80, 112 and 
160. The idea behind the matching with the ASOS model is that corrections to scaling in 
the ASOS model are similar to those of the Ising interface. Therefore it should be possible to 
obtain relaible estimates for the roughening coupling from smaller lattice sizes this way than 
from the matching with the BCSOS model. 

The total CPU requirement for all our simulations accumulates to nearly 4 years on typical 
modern workstations. The computers used include all types ranging from Pentium PCs to 
Alpha AXP-3000/400 and the parallel computer Hitachi SR2001 with 64 processors. For an 
overview of the lattice sizes employed and the CPU resources needed for the various models, 
see table [|. 

We extracted all our estimates from the matching of the two observables Di and ^3 (the 
last column in the tables). Here the convergence seems optimal. 

To obtain estimates for the roughening couplings and the matching b m we employed the 
following 

RULE: Start with the largest block lattice size, i.e., 1 = 8 (the statistical errors are the smallest 
here). As a first estimate E\ take the value for the largest lattice size L available. Then check 
whether the estimate is 2o -compatible with the results (also for 1 = 8) for the next two smaller 
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model 


lattice sizes 


CPU 


BCSOS 


16,24,32,40,48,56,64,80,96,112,128,160,192,224,256 


200 d 


XY 


32,48,64,96,192 


70 d 


ASOS 


32,48,64,96,128,192,256,384,512 


440 d 


DG 


16,24,32,48,64 


15 d 


Ising 


32,48,64,96,128,192 


650 d 



Table 1: The lattice sizes L employed in our simulations of the various models, 
together with the CPU resources needed on an "average modern workstation" . 



model 


Kr 


from I, L 




from I, L 


XY 


1.1199(1) 


4,192 


0.93(1) 


4,192 


ASOS 


0.80608(2) 


8,512 


2.78(3) 


4,512 


DG 


0.6653(2) 


2,64 


0.32(1) 


2,64 


Ising, a 


0.40759(2) 


2,192 


3.20(4) 


2,192 


Ising,b 


0.40758(1) 


8,192 


3.21(3) 


8,192 



Table 2: Our results for the roughening couplings and the matching b m , together 
with the I, L values that were used (see RULE). In the Ising model case, the index 
"a" refers to the matching with the BCSOS model, the index "6" refers to the 
matching with the ASOS model. 



lattice sizes. 2a -compatibility of two estimates mi, m<i with statistical errors e±, e<i here means 
that | mi — 77i2 1 < 2 \e\ + e^] 1 / 2 . Then also check the 2a- consistency of E\ with the estimates 
for the three largest available L-values for I = 4 and 1 = 2. If the estimates are consistent, take 
E\ as the final estimate. Otherwise restart the whole procedure with I = 4, i.e. take as the first 
estimate E\ the value from the largest L and I = 4. If there is again a failure, restart at I = 2. 

Given that the corrections die out like L~ 2 our R ULE ensures that systematic errors in the 
determination of the roughening coupling are smaller than the statistical errors quoted. 

We invite the careful reader to go through this procedure in the case of the XY-model 



(tables [l3| and 14). Our final estimates for the critical couplings and the matching b m , together 
with the values of I and L where the decision procedure stopped, are given in table ^. 

In order to determine the non-universal constants A and C we need estimates for the ratios 
of slopes R defined in eq. (p7l). These ratios for the different observables and for the different 



block/lattice sizes are presented (for the XY model as an example) in table 15. From this 
table, and from the corresponding tables for the other SOS models and the Ising model, we 
extracted by applying again our RULE a final estimate for the ratio of slopes R. Our estimates 
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model 


A 


C 


XY 


0.233(3) 


1.776(4) 


ASOS 


0.695(8) 


1.099(4) 


DG 


0.080(3) 


2.438(6 


Ising,a 


0.80(1) 


1.03(2) 


Ising,b 


0.80(1) 


1.01(1) 



Table 3: Our results for the non-universal constants A and C. In the Ising 
model case, the index "a" refers to the matching with the BCSOS model, the 
index "6" refers to the matching with the ASOS model. 



Authors 


year 




A 


C 


Swendsen p0[ 


1977 


0.77(6) 






Shugard et al. [31] 


1978 


0.68 






Janke and Nather 32] 
dito, fitl 
dito, fit2 


1991 


0.665(5) 

0.6657(3) 

0.6595(3) 


0.1204(18) 
0.0287(7) 


2.370(11) 
2.812(14) 


Evertz et al. [33] 


1993 


0.662(3) 






Hasenbusch et al. 0, || 


1992/94 


0.6645(6) 


0.078(5) 


2.44(3) 


Hasenbusch and Pinn, this work 


1996 


0.6653(2) 


0.080(3) 


2.438(6) 



Table 4: Comparison of our results for the DG model with previous estimates. 



for the non-universal constants A and C are given in table I 



5 Comparison with Previous Studies 

In this section we present a comparison of our present results with some previous estimates on 
the critical couplings Kr and non-universal parameters A and C. 

Let us start with the DG model. See table ||] for two estimates from the seventies and 
some more modern results that can be compared with the present estimates. A comparison 
of the findings in [32] and 33] with the estimates of || was presented in ref. ||. We would 
like to comment just at the apparent la-incompatibility of the present estimate for Kj^ G with 
that of |J. A closer look at our data reveals that this is most likely a statistical fluctuation: 
Discarding the L = 48 and L = 64 lattices from the analysis does not move our present estimate 
towards the result in @, which was obtained with the same method and with lattices of size 
up to L = 32. 

We now turn to the XY model. A table of previous estimates in comparison with previously 
published results is given in table H We find our present estimates consistent with our previous 
results in ||. Most of the results of the other authors, also the estimates from series analysis 
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Authors 


year 


^R 


A 


C 


Baillie and Gupta 34] 


1991 


1.1218 


0.2129 


1.7258 


Baillie and Gupta 35] 


1992 


1.119(6) 






Biferale Q36] 


1989 


1.112(2) 




1.74(20) 


Hasenbusch et al. 6| 


1992/94 


1.1197(5) 


0.223(13) 


1.78(2) 


Olsson |7 




1994 


1.12082(25) 




1.585(9) 


Olsson ]38 


1 


1995 


1.12091(13) 




1.59(2) 


Schultka and Manousakis |39J 


1994 


1.12082(16) 




1.800(2) 


Campostrini et al. |4Q| 


1996 


1.1166(4) 






Hasenbusch and Pinn, this work 


1996 


1.1199(1) 


0.233(3) 


1.776(4) 


Table 5: Comparison of our results for the XY model with previous estimates. 


Authors 


year 


tsASOS 
^R 


A 


C 


Shugard et al. g 


1978 


0.81 






Adler gl| 


1987 


0.787(24) 






Hasenbusch et al. |?], || 


1992/94 


0.8061(3) 


0.70(8) 


1.14(2) 


Hasenbusch and Pinn, this work 


1996 


0.80608(2) 


0.695(8) 


1.099(4) 



Table 6: Comparison of our results for the ASOS model with previous estimates. 



by Campostrini et al. |4C] are inconsistent with the present estimate. We conclude that in all 
these cases the systematic errors due to corrections to scaling are underestimated. 

In the case of the ASOS model, we only compare with our previous estimate || and with 
an estimate by Adler from a ninth-order low temperatur series. The series estimate has a quite 
large error, but is consistent with our result. 

We conclude this section with a comparison of the Ising interface estimates. Here we find 
that all the cited estimates of the roughening couplings are consistent with each other. Note, 
however, the large errors in the estimates that were obtained with techniques other than the 
matching method. The estimate of Mon et al. [42] for A seems to be the result of a wrong 
method. 



6 Conclusions and Outlook 

By increasing the statistics by a factor of about 100 and by also using larger lattices compared 
to ref. @, we obtained the most accurate estimates for the roughening couplings of the Ising 
interface, the ASOS model, the DG model and the 2D XY model published so far. In contrast 
to other methods systematical errors are under control. The matching procedure converges like 
L~ 2 while other methods that rely on analytic results derived from KT theory are plagued by 
corrections logarithmic in the lattice size. In addition to the precise numbers for the roughening 
coupling and other non-universal constants the matching provides an unambigous confirmation 
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Authors 


year 




A 


C 


Weeks et al. JlOj 


1973 


0.39 






Biirkner and Stauffer 11] 


1983 


0.396(22) 






Adler |tj 


1987 


0.404(12) 






Mon et al. |g] 


1988 


0.410(16) 


9.8(2.0) 


1.36(6) 


Mon et al. |3j] 


1990 


0.409(4) 






Hasenbusch ]?]] 


1992 


0.4074(3) 






Hasenbusch et al. ^ 


1996 


0.40754(5) 


0.74(2) 


1.03(2) 


Hasenbusch and Pinn, this work, a 
Hasenbusch and Pinn, this work, b 


1996 
1996 


0.40759(2) 
0.40758(1) 


0.80(1) 
0.80(1) 


1.03(2) 
1.01(1) 



Table 7: Comparison of our results for the Ising model with previous estimates. 
The index "a" referst to matching with the BCSOS model, whereas the "b" 
means matching with the ASOS model. 



of the KT nature of the phase transition of the models considered. It is interesting to compute 
the observables used for the matching method for the Sine Gordon model in perturbation 
theory. This will allow to rederive the KT flow equations from finite size scaling. Furthermore, 
it will provide quantitative information about the RG flow, in particular about the critical 
trajectory, which can be compared with the numerical results given in the present paper. 

Acknowledgements 

A large part of the simulations have been performed on the SR2001 of Hitachi Europe, Ltd. The 
remaining part of the simulations was done on workstations of the Institut fur Numerische und 
Instrumentelle Mathematik der Universitat Miinster and of the DAMTP, Cambridge University. 
This work was supported by the Leverhulme Trust under grant 16634- AOZ-R8 and by PPARC. 

References 

[1] To give only a few references: 

J. M. Kosterlitz and D. J. Thouless, J. Phys. C 6 (1973) 1181; J. M. Kosterlitz, J. Phys. 
C 7 (1974) 1046; S.T. Chui and J.D. Weeks, Phys. Rev. B 14 (1976) 4978; J.V. Jose, L.P 
Kadanoff, S. Kirkpatrick and D.R. Nelson, Phys. Rev. B 16 (1977) 1217; T. Ohta and 
K. Kawasaki, Prog. Theor. Phys. 60 (1978) 365 ; D.J. Amit, Y.Y. Goldschmidt and G. 
Grinstein, J. Phys. A 13 (1980) 585. 

[2] H. van Beijeren, Phys. Rev. Lett. 38 (1977) 993. 

[3] E.H. Lieb, Phys. Rev. 162 (1967) 162. 



13 



[4] E.H. Lieb and F.Y. Wu, in: 'Phase Transitions and Critical Phenomena', C. Domb and 
N.S. Green, eds., Vol. 1, Academic, 1972. 



[5] 
[6] 
[7] 
[8] 
[9 

[10 

[11 
[12 

[13 
[14 
[15 
[16 



[17; 

[18 

[19 
[20 
[21 
[22 
[23 



R.J. Baxter, 'Exactly Solved Models in Statistical Mechanics', Academic Press, 1982. 
M. Hasenbusch, M. Marcu, K. Pinn, Physica A 208 (1994) 124. 
M. Hasenbusch, PhD Thesis, Universitat Kaiserslautern, 1992. 



M. Hasenbusch, S. Meyer, and M. Piitz, |cond-mat/9601011 . 

H.W.J. Blote, E. Luijten, and J.R. Heringa, cond-mat/9509016, to be published in J. 
Phys. A. 

J.D. Weeks, G.H. Gilmer and H.J. Leamy, Phys. Rev. Lett. 31 (1973) 549. 

E. Burkner and D. Stauffer, Z. Phys. B 53 (1983) 241. 
M. Hasenbusch and K. Pinn, Physica A 192 (1993) 342. 

D. B. Abraham, in: 'Phase Transitions and Critical Phenomena', Vol. 10, C. Domb and 
J.L. Lebowitz, eds., Academic, 1986. 

H. van Beijeren and I. Nolden, in: 'Structure and Dynamics of Surfaces II', W. Schommers 
and P. van Blanckenhagen, eds., Topics in Current Physics Vol. 43, Springer, 1987. 

G. Forgacs, R. Lipowsky, and Th. M. Nieuwenhuizen, in: 'Phase Transitions and Critical 
Phenomena', Vol. 14, C. Domb and J.L. Lebowitz, eds., Academic, 1991. 

F. Gallet, S. Balibar, and E. Rolley, J. Physique 48 (1987) 353; F. Gallet, S. Balibar, and 

E. Rolley, J. Physique 48 (1987) 369; I.K. Robinson, E. Vlieg, H. Hornis, and E.H. Conrad, 
Phys. Rev. Lett. 67 (1991) 1890; H. Hakkinen, J. Merikoski, M. Manninen, J. Timonen, 
and K. Kaski, Phys. Rev. Lett. 70 (1993) 2451. 

R. Savit, Rev. Mod. Phys. 52 (1980) 453, and references therein. 

H. Yamamoto, Prog. Theor. Phys. 61 (1979) 363; S. Samuel, Phys. Rev. D 18 (1978) 1916; 
R.H. Swendsen, Phys. Rev. B 17 (1978) 3710; B. Nienhuis, J. Stat. Phys. 34 (1984) 731. 

I. Nolden, 'Equilibrium Crystal Shapes', PhD thesis, Utrecht 1990. 

C. Itzykson, J.M. Drouffe, 'Statistical Field Theory', Cambridge University Press 1989. 
M.P Nightingale, Physica A 83 (1976) 561. 
K. Binder, Z. Phys. B 43 (1981) 119. 
L.P. Kadanoff, Physics 2 (1966) 263. 



14 



[24] K.G. Wilson, Physica 73 (1974) 119; K.G. Wilson and J. Kogut, Phys. Rep. C 12 (1974) 
75, and references therein. 

[25] H.G. Evertz, M. Marcu, G. Lana, Phys. Rev. Lett. 70 (1993) 875. 

[26] M. Creutz, Phys. Rev. Lett. 50 (1983) 1411. 

[27] H.G. Evertz, M. Hasenbusch, M. Marcu, K. Pinn, and S. Solomon, 
Phys. Lett. B 254 (1991) 185. 

[28] M. Hasenbusch and S. Meyer, Phys. Rev. Lett. 66 (1991) 530. 

[29] M. Piitz, Diploma thesis, Kaiserslautern 1995. 

[30] R.H. Swendsen, Phys. Rev. B 15 (1977) 5421. 

[31] W.J. Shugard, J.D. Weeks, and G.H. Gilmer, Phys. Rev. Lett. 41, 1399 (1978). 

[32] W. Janke and K. Nather, Phys. Lett. A 157, 11 (1991). 

[33] H.G. Evertz, M. Hasenbusch, M. Marcu, K. Pinn, Physica A 199 (1993) 31. 

[34] C.F. Baillie and R. Gupta, Nucl. Phys. B (Proc. Suppl.) 20 (1991) 669. 

[35] C.F. Baillie and R. Gupta, Phys. Rev. B 45 (1992) 2883. 

[36] L. Biferale and R. Petronzio, Nucl. Phys. B 328 (1989) 677. 

[37] P. Olsson, Phys. Rev. Lett. 73 (1994) 3339. 

[38] P. Olsson, Phys. Rev. B 52 (1995) 4526. 

[39] N. Schultka and E. Manousakis, Phys. Rev. B 49 (1994) 12071. 

[40] M. Campostrini, A. Pelisetto, P. Rossi and E. Vicari, [hep-lat/9603002 . 

[41] J. Adler, Phys. Rev. B 36 (1987) 2473. 

[42] K.K. Mon, S. Wansleben, D.P Landau and K. Binder, Phys. Rev. Lett. 60 (1988) 708. 
[43] K.K. Mon, D.P. Landau and D. Stauffer, Phys. Rev. B 42 (1990) 545. 



15 



List of Tables 



1 Lattices sizes and CPU requirements 

2 Results for the roughening couplings and the matching b. 

3 Results for the non-universal constants A and C ... 

4 Comparison of our DG results with previous estimates 

5 Comparison of our XY results with previous estimates 

6 Comparison of our ASOS results with previous estimates 

7 Comparison of our Ising results with previous estimates 

8 Exact results for A\ ^ and 

9 BCSOS results for the critical coupling, I = 1, 2 

10 BCSQS results for the critical coupling, I = 4, 8 

11 Slope of BCSOS observables for the critical coupling, / = 1, 2 

12 Slope of BCSOS observables for the critical coupling, / = 4, 8 

13 XY results for the matching factor 

14 XY results for the roughening coupling 

15 Ratio of slopes for the XY model 



16 



T 

LJ 


^1,2 


,(o) 

^2,2 


4 (o) 


4(0) 


^1,8 


AO) 
^2,8 


lb 





1231474 





1718750 


U.243918U 


0.326U9U5 





3177271 


U.4zo014o 


94 





1205651 





1689815 


D 9^41 4fifi 







2815115 


U.OOO^iUUU 


32 





1196545 





1679687 


0.2306619 


0.3119780 





2682336 


0.3691650 


40 





1192317 





1675000 


0.2290354 


0.3102652 





2619443 


0.3624936 


48 





1190016 





1672454 


0.2281480 


0.3093327 





2584828 


0.3588444 


56 





1188628 





1670918 


0.2276114 


0.3087697 





2563784 


0.3566345 


64 





1187726 





1669922 


0.2272625 


0.3084039 





2550051 


0.3551961 


80 





1186664 





1668750 


0.2268516 


0.3079734 





2533820 


0.3535002 


96 





1186088 





1668113 


0.2266280 


0.3077394 





2524963 


0.3525768 


112 





1185740 





1667730 


0.2264931 


0.3075982 





2519608 


0.3520192 


128 





1185514 





1667480 


0.2264055 


0.3075066 





2516126 


0.3516570 


160 





1185248 





1667187 


0.2263024 


0.3073988 





2512025 


0.3512307 


192 





1185104 





1667028 


0.2262464 


0.3073402 





2509794 


0.3509989 


224 





1185017 





1666932 


0.2262126 


0.3073049 





2508448 


0.3508592 


256 





1184960 





1666870 


0.2261907 


0.3072820 





2507573 


0.3507684 


384 





1184858 





1666757 


0.2261509 


0.3072403 





2505986 


0.3506036 


512 





1184822 





1666718 


0.2261370 


0.3072258 





2505429 


0.3505459 


oo 





118478 





166667 


0.226119 


0.307207 





250471 


0.350472 



Table 8: Exact results for A\ and A 2 as functions of the size of the funda- 
mental lattice (L) and the size of the blocked lattice (7). The last row contains 
values extrapolated to L = 00. 
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L 


I 




A 2 


A 3 


A 4 


16 


1 






0.26511 


;30) 


0.06737 


'26) 


24 


1 






0.24004 


;3i) 


0.05437 


[25) 


32 


1 






0.22572 


;32) 


0.04772 


[25) 


40 


1 






0.21526 


;3i) 


0.04341 


[24) 


48 


1 






0.20778 


;30) 


0.03959 


[24) 


56 


1 






0.20157 


;3i) 


0.03704 


[24) 


64 


1 






0.19656 


;33) 


0.03571 


[24) 


80 


1 






0.18884 


;33) 


0.03252 


[25) 


96 


1 






0.18209 


;34) 


0.03032 


[25) 


112 


1 






0.17797 


;35) 


0.02880 


[25) 


128 


1 






0.17399 


;33) 


0.02750 


[23) 


160 


1 






0.16704 


;37) 


0.02529 


[26) 


192 


1 






0.16341 


10) 


0.02419 


[26) 


224 


1 






0.15935 


ii) 


0.02265 


[26) 


256 


1 






0.15616 


;42) 


0.02188 


[26) 


16 


2 


0.085740(27) 


0.117435(45) 


0.22382 


18) 


0.05597 


13) 


24 


2 


0.083135(26) 


0.115240(43) 


0.20229 


18) 


0.04436 


12) 


32 


2 


0.082161(25) 


0.114434(42) 


0.18981 


19) 


0.03851 


12) 


40 


2 


0.081661(24) 


0.113976(41) 


0.18079 


19) 


0.03444 


12) 


48 


2 


0.081351(24) 


0.113660(41) 


0.17409 


18) 


0.03174 


12) 


56 


2 


0.081129(24) 


0.113497(41) 


0.16888 


19) 


0.02975 


12) 


64 


2 


0.080932(24) 


0.113249(41) 


0.16449 


'20) 


0.02812 


12) 


80 


2 


0.080719(24) 


0.113064(41) 


0.15767 


[20) 


0.02567 


12) 


96 


2 


0.080541(24) 


0.112855(41) 


0.15218 


[21) 


0.02401 


12) 


112 


2 


0.080422(24) 


0.112743(41) 


0.14834 


[21) 


0.02275 


12) 


128 


2 


0.080306(22) 


0.112617(38) 


0.14485 


[20) 


0.02145 


I") 


160 


2 


0.080165(25) 


0.112424(42) 


0.13916 


[23) 


0.01983 


12) 


192 


2 


0.080055(25) 


0.112278(43) 


0.13567 


[24) 


0.01848 


12) 


224 


2 


0.080008(26) 


0.112258(44) 


0.13225 


[25) 


0.01763 


12) 


256 


2 


0.079953(26) 


0.112235(45) 


0.12936 


[26) 


0.01699 


13) 



Table 9: Monte Carlo results for the Ai obtained at the critical couplin; 
the BCSOS model. The block-observables Ai are defined in the text. L is 
original lattice size, and / is the size of the blocked system. 
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L 


I 


A 


1 




A 2 


A 


3 


A 


4 




16 


4 





176651 


[27) 


0.230216(40) 





221936(80) 





077948 


[65) 


24 


4 





165308 


;24) 


0.219807(37) 





194147(82) 





052961 


'63) 


32 


4 





161166 


;23) 


0.215927(35) 





179384(82) 





043583 


;ei) 


40 


4 





159116 


[22) 


0.213952(35) 





169635(82) 





038175 


[60) 


48 


4 





157915 


[22) 


0.212736(34) 





162291(81) 





034469 


[59) 


56 


4 





157067 


[22) 


0.211902(34) 





156869(82) 





031988 


[59) 


64 


4 





156463 


[21) 


0.211246(34) 





152173(85) 





029793 


[59) 


80 


4 





155714 


[21) 


0.210486(34) 
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[58) 


96 


4 
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024779 
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4 
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[21) 


0.209486(34) 
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4 
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126456(96) 





019967 


[59) 


192 
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153818 
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0.208360(35) 
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018579 
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4 





153583 


[22) 


0.208080(35) 
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017534 
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4 





153391 


[23) 
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016720 


[59) 


16 


8 





265665 


;22) 


0.325828(30) 





345213(51) 


1 
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24 


8 





216030 


16) 


0.280733(24) 





227866(37) 





105810 


[34) 


32 


8 





196646 


14) 


0.264669(22) 





205656(38) 





072017 


[32) 


40 


8 





188911 


14) 


0.257187(21) 





189140(37) 





055770 


[31) 


48 


8 





184346 


13) 


0.252912(20) 





179067(36) 





048696 


[30) 


56 


8 





181606 


13) 


0.250139(20) 





171361(37) 





043493 


[30) 


64 


8 





179710 


13) 


0.248261(20) 





165313(38) 





039956 


[30) 


80 


8 





177383 


12) 


0.245856(20) 





156139(38) 





034926 


[29) 


96 


8 





175951 


12) 


0.244337(20) 





149442(40) 





031601 


[30) 


112 


8 





174995 


12) 


0.243293(20) 





144178(39) 





029200 


[29) 


128 


8 





174309 


11) 


0.242532(18) 





139959(37) 





027326 


[27) 


160 


8 





173355 


12) 


0.241458(20) 





133345(41) 





024585 


[29) 


192 


8 





172733 


12) 


0.240736(20) 





128466(43) 





022659 


[29) 


224 


8 





172270 


12) 


0.240186(20) 





124527(44) 





021187 


[29) 


256 


8 





171898 


13) 


0.239742(20) 





121338(46) 





020028 


[29) 



Table 10: Continuation of table |9|. 
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A 3 


A 4 


16 
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1 
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2.286(57) 


192 


1 
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224 


1 






8.478(81) 


2.432(79) 
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1 






8.719 95 


2.348 91) 


16 


2 


-0.54495(72) 
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1.6789(32) 


24 


2 


-0.56652(96) 


-0.8146(16) 


4.4235(56) 


1.6717 47 


32 


2 


-0.5876(12) 


-0.8476(20) 


4.7222(73) 


1.6559(59) 


40 


2 


-0.6055(15) 


-0.8725 25 

V / 


4.9708(89) 


1.6655 74 


48 


2 


-0.6208(17) 


-0.8961(29) 


5.172(11) 


1.6574(88) 


56 


2 


-0.6331(20) 


-0.9126(34) 


5.345(13) 


1.659(10) 


64 


2 


-0.6469(23) 


-0.9316(38) 


5.484(14) 


1.642(12) 


80 


2 


-0.6644(28) 


-0.9587(47) 


5.744(18) 


1.651(14) 


96 


2 


-0.6806(33) 


-0.9807(56) 


5.907(21) 


1.620(17) 


112 


2 


-0.6822(39) 


-0.9882(67) 


6.056(25) 


1.659(20) 


128 


2 


-0.7079(40) 


-1.0227(68) 


6.189(26) 


1.639(21) 


160 


2 


-0.7199(55) 


-1.0320(94) 


6.474(36) 


1.641(29) 


192 


2 


-0.7412(64) 


-1.072(11) 


6.571(44) 


1.642(35) 


224 


2 


-0.7613(75) 


-1.103(13) 


6.734(50) 


1.655(40) 


256 


2 


-0.7748(88) 


-1.121(15) 


6.939(58) 


1.592(46) 



Table 11: Monte Carlo results for the derivatives of the Ai with respect to 
the coupling, taken for the BCSOS model at the critical coupling. The block- 
observables Ai are defined in the text. L is the original lattice size, and / is the 
size of the blocked system. 
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L 


I 


A x 


A 2 


A 3 


A 4 


16 


4 


-1.02443(81) 


-1.3944(12) 


2.9483(20) 


1.5854(17) 


24 


4 


-1.0343(10) 


-1.4308(15) 


3.2837(29) 


1.4291(24) 


32 


4 


-1.0626(12) 


-1.4790(18) 


3.5337(36) 


1.3832(30) 


40 


4 


-1.0909(14) 


-1.5208(22) 


3.7292(44) 


1.3796(37) 


48 


4 


-1.1165(17) 


-1.5573(26) 


3.8885(51) 


1.3631(44) 


56 


4 


-1.1395(19) 


-1.5901(30) 


4.0296(61) 


1.3619(51) 


64 


4 


-1.1596(21) 


-1.6176(33) 


4.1442(67) 


1.3456(58) 


80 


4 


-1.1943(26) 


-1.6676(41) 


4.3474(85) 


1.3524(72) 


96 


4 


-1.2182(31) 


-1.7022(48) 


4.491(10) 


1.3459(87) 


112 


4 


-1.2410(36) 


-1.7339(56) 


4.627(12) 


1.3395(99) 


128 


4 


-1.2668(37) 


-1.7680(58) 


4.744(12) 


1.329(10) 


160 


4 


-1.2935(50) 


-1.7997(78) 


4.942(17) 


1.361(14) 


192 


4 


-1.3372(59) 


-1.8660(92) 


5.055(20) 


1.372(17) 


224 


4 


-1.3603(68) 


-1.897(11) 


5.186(23) 


1.352(20) 


256 


4 


-1.3771(79) 


-1.921(12) 


5.351(27) 


1.330(23) 


16 


8 


-1.27254(71) 


-1.65259(97) 


2.6102(15) 


0.00000 


24 


8 


-1.09526(73) 


-1.5480(11) 


2.4655(15) 


1.4653(13) 


32 


8 


-1.09968(81) 


-1.5481(12) 


2.7517(19) 


1.4710(16) 


40 


8 


-1.10306(92) 


-1.5727(14) 


2.9006(22) 


1.3380(19) 


48 


8 


-1.1188(10) 


-1.6001(16) 


3.0496(26) 


1.3139(22) 


56 


8 


-1.1344(12) 


-1.6265(18) 


3.1721(30) 


1.2878(26) 


64 


8 


-1.1514(13) 


-1.6537(20) 


3.2715(34) 


1.2758(30) 


80 


8 


-1.1828(16) 


-1.7011(25) 


3.4528(41) 


1.2607(36) 


96 


8 


-1.2067(18) 


-1.7356(29) 


3.5904(50) 


1.2524(43) 


112 


8 


-1.2338(21) 


-1.7768(33) 


3.7300(58) 


1.2482(50) 


128 


8 


-1.2553(21) 


-1.8070(34) 


3.8241(59) 


1.2377(52) 


160 


8 


-1.2906(29) 


-1.8574(46) 


4.0041(81) 


1.2439(71) 


192 


8 


-1.3261(34) 


-1.9101(54) 


4.1256(99) 


1.2387(85) 


224 


8 


-1.3503(39) 


-1.9428(62) 


4.252(11) 


1.240(10) 


256 


8 


-1.3741(45) 


-1.9792(72) 


4.382(13) 


1.230(11) 



Table 12: Continuation of table |ll]. 



21 



L 


I 




A 2 ,As 




D 2 ,A 3 


32 


2 


0.9176(46) 


0.9381(65) 


0.9007(55) 


0.9273(74) 


48 


2 


0.9239(62) 


0.9307(78) 


0.9148(69) 


0.9235(86) 


64 


2 


0.9206(80) 


0.9275(92) 


0.9148(86) 


0.9231(97) 


96 


2 


0.9362(97) 


0.9386(117) 


0.9335(101) 


0.9365(121) 


128 


2 


0.9376(107) 


0.9356(118) 


0.9362(109) 


0.9345(120) 


192 


2 


0.9275(136) 


0.9331(170) 


0.9266(138) 


0.9324(172) 


32 


4 


0.8896(14) 


0.9035(17) 


0.8601(18) 


0.8820(20) 


48 


4 


0.9080(20) 


0.9162(23) 


0.8923(24) 


0.9045(26) 


64 


4 


0.9172(27) 


0.9231(31) 


0.9082(29) 


0.9165(33) 


96 


4 


0.9263(36) 


0.9300(38) 


0.9217(38) 


0.9266(40) 


128 


4 


0.9364(40) 


0.9363(43) 


0.9340(41) 


0.9345(45) 


192 


4 


0.9287(50) 


0.9271(52) 


0.9273(51) 


0.9260(52) 


32 


8 


0.8421(3) 


0.8622(4) 


0.7731(6) 


0.8149(6) 


48 


8 


0.8760(6) 


0.8903(7) 


0.8329(8) 


0.8610(9) 


64 


8 


0.8890(8) 


0.9009(10) 


0.8629(10) 


0.8825(11) 


96 


8 


0.9118(11) 


0.9199(13) 


0.8988(13) 


0.9109(14) 


128 


8 


0.9185(14) 


0.9245(16) 


0.9109(16) 


0.9191(17) 


192 


8 


0.9267(19) 


0.9297(21) 


0.9227(20) 


0.9269(22) 



Table 13: XY results for the matching factor obtained in the way described 
after eq. (§§). 
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L 


I 


Ai,A 3 


-42,-43 




D 2 Al 


32 


2 


1.119982(94) 


1.119545(117) 


1.120352(105) 


1.119773(133) 


48 


2 


1.119933(90) 


1.119822(120) 


1.120084(98) 


1.119940(130) 


64 


2 


1.119757(95) 


1.119661(114) 


1.119839(100) 


1.119722(119) 


96 


2 


1.119707(86) 


1.119682(101) 


1.119735(89) 


1.119703(104) 


128 


2 


1.119831(92) 


1.119852(102) 


1.119845(94) 


1.119863(104) 


192 


2 


1.119789(92) 


1.119742(114) 


1.119797(93) 


1.119748(115) 


32 


4 


1.120176(48) 


1.119726(53) 


1.121162(52) 


1.120428(59) 


48 


4 


1.120034(46) 


1.119846(51) 


1.120403(50) 


1.120116(55) 


64 


4 


1.119868(45) 


1.119754(52) 


1.120044(47) 


1.119883(55) 


96 


4 


1.119791(41) 


1.119735(49) 


1.119860(42) 


1.119785(51) 


128 


4 


1.119798(41) 


1.119800(48) 


1.119829(42) 


1.119823(48) 


192 


4 


1.119866(42) 


1.119884(47) 


1.119881(43) 


1.119896(48) 


32 


8 


1.120078(32) 


1.118559(33) 


1.125145(34) 


1.122223(39) 


48 


8 


1.120069(29) 


1.119442(29) 


1.121941(28) 


1.120741(33) 


64 


8 


1.120095(26) 


1.119701(28) 


1.120979(28) 


1.120310(30) 


96 


8 


1.119951(26) 


1.119754(26) 


1.120271(26) 


1.119973(27) 


128 


8 


1.119906(24) 


1.119789(26) 


1.120057(25) 


1.119894(27) 


192 


8 


1.119888(22) 


1.119841(24) 


1.119949(23) 


1.119885(24) 



Table 14: XY results for the roughening coupling obtained in the way described 
after eq. (§§). 
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L 


I 




D 2 


Ax 


A 2 


16 


1 






-0.4167(14) 


-0.4137(17) 


24 


1 






-0.4206(14) 


-0.4195(25) 


32 


1 






-0.4249(17) 


-0.4227(33) 


48 


1 






-0.4252(19) 


-0.4267(48) 


64 


1 






-0.4255(24) 


-0.4270(67) 


96 


1 






-0.4257(32) 


-0.4410(99) 


128 


1 






-0.4267(41) 


-0.439(13) 


192 


1 






-0.4280(59) 


-0.413(19) 


16 


2 


-0.4271(11) 


-0.4302(11) 


-0.4154(13) 


-0.4046(12) 


24 


2 


-0.4299(13) 


-0.4314(13) 


-0.4194(13) 


-0.4112(17) 


32 


2 


-0.4298(15) 


-0.4304(14) 


-0.4236(14) 


-0.4167(23) 


48 


2 


-0.4288(19) 


-0.4285(19) 


-0.4244(17) 


-0.4195(34) 


64 


2 


-0.4281(23) 


-0.4285(24) 


-0.4238(20) 


-0.4271(46) 


96 


2 


-0.4355(33) 


-0.4365(33) 


-0.4263(25) 


-0.4326(70) 


128 


2 


-0.4250(41) 


-0.4252(42) 


-0.4281(32) 


-0.4140(92) 


192 


2 


-0.4282(60) 


-0.4277(62) 


-0.4276(46) 


-0.412(14) 


16 


4 


-0.41961(73) 


-0.42505(77) 


-0.4070(13) 


-0.3657(12) 


24 


4 


-0.42592(89) 


-0.42839(90) 


-0.4140(13) 


-0.3962(11) 


32 


4 


-0.42965(98) 


-0.43038(95) 


-0.4182(13) 


-0.4061(14) 


48 


4 


-0.4290(12) 


-0.4291(11) 


-0.4218(13) 


-0.4148(21) 


64 


4 


-0.4294(13) 


-0.4296(13) 


-0.4213(15) 


-0.4187(28) 


96 


4 


-0.4324(17) 


-0.4323(17) 


-0.4250(17) 


-0.4236(43) 


128 


4 


-0.4296(21) 


-0.4290(20) 


-0.4264(21) 


-0.4177(57) 


192 


4 


-0.4295(30) 


-0.4300(29) 


-0.4281(29) 


-0.4050(85) 


24 


8 


-0.42239(92) 


-0.42300(69) 


-0.4120(18) 


-0.38303(53) 


32 


8 


-0.42373(70) 


-0.42791(81) 


-0.4087(11) 


-0.3694(18) 


48 


8 


-0.42871(87) 


-0.43010(86) 


-0.4147(13) 


-0.3952(12) 


64 


8 


-0.42960(99) 


-0.43030(95) 


-0.4177(13) 


-0.4035(15) 


96 


8 


-0.4320(12) 


-0.4328(12) 


-0.4204(15) 


-0.4159(23) 


128 


8 


-0.4313(14) 


-0.4311(13) 


-0.4228(15) 


-0.4192(31) 


192 


8 


-0.4311(18) 


-0.4307(17) 


-0.4267(19) 


-0.4175(48) 



Table 15: Monte Carlo estimates for the ratio of slopes for the XY and BCSOS 
model block observables, taken at the critical couplings, c.f. the definition in 
eq. ©. 
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